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A novel spin-spin coupling mechanism that occurs during the transport of spin-polarized minority 
electrons in semiconductors is described. Unlike the Coulomb spin drag, this coupling arises from 
the ambipolar electric field which is created by the differential movement of the photoelectrons and 
the photoholes. Like the Coulomb spin drag, it is a pure spin coupling that does not affect charge 
diffusion. Experimentally, the coupling is studied in GaAs using polarized microluminescence. 
The coupling manifests itself as an excitation power dependent reduction in the spin polarization 
at the excitation spot without any change of the spatially averaged spin polarization. 


Modification of diffusion to include ambipolar effects in 
charged heterogeneous media is a topic of interest in sev¬ 
eral fields, including the study of astrophysical objects [l| 
and plasmas [2|, as well as in semiconductors The 

case of semiconductors is of interest both for applications 
[il and because electron gases in semiconductors can be 
spin polarized. The question of the effect of ambipolar 
coupling on spin polarized carrier diffusion is yet to be 
addressed 0, but is likely to be of importance for any 
future bipolar semiconductor spintronic device. A vari¬ 
ety of spin-spin and/or spin-charge coupling phenomena 
have been revealed in semiconductor^ for example those 
due to the spin-orbit interaction [3, [H , those due to the 
Pauli principle @,[ 13 , as well as the Coulomb spin drag 
[ill , fl^ . In the latter case it was shown that a coupling 
between the + and — spins results in a spin diffusion con¬ 
stant whose magnitude is smaller than the charge diffu¬ 
sion constant. 

Here we describe and study a novel spin-spin coupling 
mechanism of ambipolar origin which yields a spin diffu¬ 
sion constant whose magnitude is larger than the charge 
diffusion constant. This coupling occurs in the presence 
of a spatially inhomogeneous gas of spin-polarized photo¬ 
electrons and of unpolarized, slower diffusing holes. The 
differential diffusion of + spin electrons and holes cre¬ 
ates an internal electric field which acts on both + and 

— spins thereby coupling them. The same is true for the 

— spin electrons. A full description of the coupled dif¬ 
fusion equations is given and the effect is experimentally 
observed in GaAs. 

The sample is a 3 /im thick, Be-doped {Na = 1.5 x 
10^^ cm“^) GaAs film covered on both sides by passivat¬ 
ing GalnP layers which not only reduce the surface re¬ 
combination velocity, but confine the photocarriers to the 
active layer. The sample is studied using a microlumines¬ 
cence technique described elsewhere [IJ that, as shown 
in Fig. [U creates a spatially inhomogeneous population 
of spin polarized electrons and of unpolarized holes. This 
is a pre-requisite for the observation of ambipolar cou¬ 
pling phenomena. The photoexcitation is achieved us¬ 
ing a tightly-focused circularly-polarized GW pump (1/e 
half width oi w = 0.6 /im, energy 1.59 eV) so that, at 
the chosen value of Na^ ambipolar coupling becomes im¬ 


portant for experimentally accessible pump powers. All 
experiments reported here are at 300 K where other cou¬ 
pling phenomena I,[Mm are negligible. The lumines¬ 
cence intensity and polarization are monitored as a func¬ 
tion of distance, r, from the excitation spot, from which 
depth integrated profiles of the photoelectron charge den¬ 
sity n = + n_, and the photoelectron spin density 

s = — n_, can be obtained respectively [T3|. Here 

n± are the concentrations of electrons of spin ± with 
a quantization axis chosen along the direction of light 
propagation. 



FIG. 1. (a) The principle of the experimental technique 

in which the sample is photoexcited by a tightly-focused, 
circularly-polarized pump (red arrows and lines), and the re¬ 
sulting photoluminescence intensity (blue arrows and lines) 
and polarization is measured. These two quantities can be 
used to obtain depth (z) integrated charge and spin density 
profiles as a function of r, the radial distance from the exci¬ 
tation spot, (b) Room temperature charge and spin density 
images obtained from the GaAs sample at low excitation 
power (0.03 mW) where ambipolar coupling is negligible, (c) 
The angular averaged profiles obtained from the images in (b). 
The solid lines are fits using the solution to the uncoupled dif¬ 
fusion equations 0, from which one obtains the charge (L) 
and spin (Lg) diffusion lengths. 

The charge and spin density profiles at low power 
(0.03 mW) illustrate the unipolar regime. These profiles, 
shown in Fig. HTc) . ar e analyzed using the uncoupled 
diffusion equations [I^ (solid lines in the figure), from 
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which the charge diffusion length L = =10 /rm 

and the spin diffusion length Lg = \/DeTs = 0.95 /im 
are obtained. The spatially averaged spin polarization, 
defined as<^>=<5>/ <n>, does not depend on 
diffusion and is equal, in a two-dimensional picture, to 
V)^ where is the effective initial polarization 
including possible losses during thermalization or during 
diffusion along the z direction. The experimental values 
of < ^ 0.4% shown in Fig. EJa), along with the val¬ 

ues of L and imply = 0.4, slightly smaller than 
its value of 0.5 without losses 0 . The transport param¬ 
eters for the sample are then characterized by assuming 
an electron mobility of /ie = 3350 cm^/Vs for this doping 
level 00, from which the charge diffusion coefficient, 
De = 86 cm^/s, is obtained using the Einstein relations. 
Combining this with the measured value of L, a minority 
carrier lifetime r = I? jD^ = 11.6 ns is found, close to 
that measured in similarly doped GaAs [5|. Assuming 
that in the unipolar limit the spin diffuses with the same 
diffusion coefficient as the charge, the value of implies 
a spin relaxation time of T\ = 105 ps. The hole mobility 
is assumed to be = 220 cm^/Vs 0 . As will be seen 
below, these are all the parameters necessary to describe 
the spin-spin coupling observed at higher excitation pow¬ 
ers. 
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FIG. 2. The spatial profiles of the spin polarization ^ — sjn 
with increasing values of the excitation power along with the 
corresponding power dependence of the spatially-averaged po¬ 
larization < ^ >—< s > / < n > (upper inset) shown as 
obtained experimentally (a), calculated using Eq. [3 which in¬ 
cludes the spin-spin coupling (b), and calculated using Eq. [8] 
which does not include this coupling. The experimental data 
reveal a decrease in ^ with increasing excitation power at 
r = 0 by about a factor of 6 without any significant change 
in < ^ >. A comparison of the experimental and calcu¬ 
lated curves reveals that the large reduction in ^(r = 0) is a 
manifestation of the ambipolar spin-spin coupling. 


Fig. Efa) shows the spin polarization profiles for 
increasing excitation powers. At low power, one has 
^{r = 0) = 11% which is a factor 25 larger than < > 

because the effective lifetime at r = 0 is not r but a 
diffusion time {rdiff ~ j^D^. Since r^iff ~ 10 ps 

^ r, this reduces polarization losses by relaxation at 
r = 0. Upon increasing the pump power, the polariza- 
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FIG. 3. (a) The measured polarization at r = 0 along with 
the spatially averaged polarization < > as a function of 

excitation power in the GaAs sample doped at Na — 1.5 x 10^^ 
cm“^. Also shown are the polarizations at r = 0 calculated 
with (black line) and without (green line) spin-spin coupling, 
(b) The same data shown for a GaAs sample with Na — 
10^® cm“^, for which ambipolar effects are strongly reduced 
and the spin-spin coupling is absent. 


tion at r = 0 decreases by almost a factor of 6 to ^ 2% at 
the maximum accessible power (3 mW). It is important 
to note that this decrease is not due to a decrease in T\ 
since the spatially-averaged polarization does not change 
with excitation power (see upper panel of Fig. Efa)). 

In order to interpret the experimental results it is nec¬ 
essary to calculate the internal ambipolar fields, cre¬ 
ated by diffusion of electrons of ± spins. Neglecting the 
effects of Pauli blockade, thermoelectric phenomena, and 
spin Goulomb drag which is screened by the majority 
holes [T^ , the conservation equations for + spins and for 
holes are: 


^+-n+/T-(n+-n_)/(2Ti)-kV-[cr+E/g^ + L)eVn+] = 0 

( 1 ) 

(^+ + 9-) - + V • [-ahE/qN Dh^Sp] = 0, (2) 


where the conservation equation for — spins is obtained 
by exchanging + and — in Eq. [TJ The generation rate 
g± of ± spins is strongly peaked at r = 0 since a tightly 
focused light excitation is used. Here, 6p is the photohole 
concentration E is the internal ambipolar electric field, q 
is the absolute value of the electron charge and is the 
hole diffusion constant. The spin conductivities are given 
by CF± = qpe'^± and the electron and hole conductivi¬ 
ties are respectively, (Jq = qpe"^ and (Jh = qjiih{E^ + Sp), 
where is the concentration of charged acceptors. Gal- 
culation of the electric field will assume a 2-dimensional 
picture (to be justified below) in which the divergences of 
Eq. [1] and Eq. [2] are approximated by derivatives in the 
sample plane. Gombination of Eq. [T]for + spins with the 
same equation for — spins and with Eq. [2] for holes shows 
that E can be written as the sum of three contributions: 


- DhiVSp - Vn) 
= q - 




( 3 ) 
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The first term is caused by a possible disruption of lo¬ 
cal charge neutrality {Sp ^ n) and will be assumed, as 
verified below, to be negligible. The contribution Ej^ is 
given by: 

E+ = q \ W n+. (4) 

This field is proportional to the difference in the diffusiv- 
ities of electrons and holes and is identified as the am- 
bipolar field generated by diffusion of - 1 - spins. A similar 
expression is obtained for E- that is associated with — 
spins. 

The current J+ of spins + is then the sum of the dif¬ 
fusion current and of the drift current in Ej^ + E _. The 
spin-spin coupling is explicit since for example can be 
decomposed into two components, J++ and J+-, propor¬ 
tional to Vn+ and to Vn_, respectively. This current, 
together with J_, is given by 


4 


= q 


[dp 




(5) 


where (erg + ah)D'^~^ = a^Dh + {(Jh Ecr-)De and the non 
diagonal element, given by {ae-\-cFh)D^~ = a-^{Dh— De) 
is zero in the unipolar case where De = Dh- The other 
matrix elements are obtained by exchanging + and —. 
The form of the diffusion matrix is similar to that de¬ 
scribing the Coulomb spin drag , with the notable dif¬ 
ference that the nondiagonal elements are negative since 
an outward diffusion of - 1 - spins generates an inward fiux 
of — spins. 

It is pointed out that, for a nonzero spin polarization, 
the two nondiagonal elements of the diffusion matrix are 
not equal. However, the corresponding diffusion currents, 

P_=q{Dh-De)^^^. ( 6 ) 

O'e + 


and J _^ (obtained by exchanging + and —) are equal 

to first order. Indeed, for nondegenerate electrons, the 
diffusion length is the same for the two types of spins, so 

that Vn+/Vn_ ^ n+/n_. The currents _and J_+ 

describe the flow of comparable numbers of - 1 - and — spins 
per unit time towards the excitation spot which therefore 
reduces the polarization at r = 0. This is indeed what is 
observed experimentally. 

Further insight into the nature of the ambipolar spin- 
spin coupling is gained by studying the diffusion equa¬ 
tions for n and 5 , obtained from Eq. [Uby replacing E by 
its value defined in Eq. [3l It is immediately clear that 
the coupling described by the off-diagonal terms of Eq. [5] 
is a pure spin effect, since the charge diffusion equation 
becomes: 


( 5 ++ 5 -) - V'T + V • [DaVn] = 0, (7) 

where the unipolar diffusion constant is replaced by the 
usual ambipolar diffusion one Q, Da defined by (de + 


(Th)Da = {cTeDh + (JhDe)^ without any coupling to the 
electronic spins. On the other hand spin-spin coupling 
modifies the spin conservation equation which becomes: 

{g+ - 5-) - + V ■ [{D, -D',ps + D', ^Vn] = 0 , ( 8 ) 


where 


D'^ = Dt-+D- 


-+ 


ae{Dh - De) 


(9) 


If spin relaxation is negligible the electronic spin polar¬ 
ization is spatially constant (implying ^Vn = Vs) and 
the divergence term in Eq. [8] reduces to V • [DaVs]. Spin 
then diffuses in the same way as charge. In the oppo¬ 
site case where ^ decreases with distance as is generally 
found for local light excitation the divergence term 
of Eq. [8] is of the form V • [{Da — D'^){Vs + 6^Vn)], 
where S = D'^jD^ is close to -1 at high power. Two im¬ 
portant conclusions are to be drawn from this analysis. 
Eirstly, the spin diffusion constant Da — D'^ is now larger 
than the charge diffusion constant which is in direct con¬ 
trast with Coulomb spin drag. Secondly, spin diffusion 
now depends on charge due to the 6^Vn term. This 
term has the same form as that induced by Coulomb 
spin drag, or by diffusion of degenerate spins, with the 
notable difference that in these two cases > 0 00. 

The agreement between the above model and the ex¬ 
perimental results is now verified quantitatively using a 
numerical resolution of Eq. 0 and Eq. [8] with the sam¬ 
ple parameter values determined above. As shown in Eig. 
ma), in a low power unipolar regime one finds a polariza¬ 
tion at r = 0 of 11.5%, very close to the measured value. 
Using D'^ = 0 (no spin-spin coupling), one finds a slight 
decrease of polarization at r = 0 with increasing power 
which is not sufficient to explain the experimental obser¬ 
vations. This decrease is mostly due to the decrease of 
the ambipolar diffusion constant Da which results in an 
increase of the diffusion time Tdiff. On the other hand, 
the inclusion of the spin-spin coupling term, as shown 
in Eig. El^b) accounts very well for the experimental re¬ 
sults. It is Anally verified that Pauli blockade 0 does 
not play a role here; for the maximum power the depth 
averaged value of n(r = 0) is calculated to be of the or¬ 
der of n = 1.6 X 10^^ cm“^, which is smaller than the 
effective density of states of the conduction band at 300 
K. 

The comparison between experiment and theory is 
summarized in Eig. Ha) which shows excellent agree¬ 
ment with the experimentally measured power depen¬ 
dence of the spin polarization at r = 0, and the poor 
one obtained if D'^ = 0. The ambipolar nature of the ex¬ 
perimental results is further confirmed by measurements 
at 300 K using another sample with an increased accep¬ 
tor doping of Na = 10^^ cm“^. At this doping density 
ambipolar coupling is strongly reduced because of the in¬ 
creased majority hole conductivity, so that Da ~ Dq and 
—D'^ Dq. Indeed, as shown in Eig. [31(b), the polariza¬ 
tion at the excitation spot does not exhibit any decrease 
as a function of excitation power. 
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FIG. 4. (a) Spatial dependence of the relative difference be¬ 
tween the photoelectron and photohole density for increasing 
values of the excitation power, (b) Internal ambipolar elec¬ 
tric field. In spite of the small relative difference between the 
local electron and hole charge densities, this field can be quite 
large near r = 0 at high excitation power, (c) The electric 
field obtained from the approximate expression Eq. |3]is close 
to the numerically calculated result in (b). 


The two main approximations - local charge neu¬ 
trality and the two-dimensional nature of the diffusive 
transport - are now justified. The hypothesis of local 
charge neutrality is verified by removing the approxima¬ 
tion n = 6p and by performing a numerical resolution 
of Eq. [1] and Eq. [2l together with Poisson’s equation 
V ' E = {q/eeo){6p — n), where cq is the vacuum permit¬ 
tivity and e is the dielectric constant of GaAs. As shown 
in Fig. m^a), it is found that the relative photoinduced 
electric charge | (n — Sp) / {n -\- Sp) \ is always smaller than 


10“^. As shown in Fig. SJ^b), the internal electric field is 
peaked near r = 1 pm and can be as large 100 V/cm at 
high power. This corresponds to a drift length of several 
tens of pm i.e. larger than the charge diffusion length 
and implies that near r = 0 charge and spin drift in the 
internal electric field prevails over diffusion. This is the 
origin of the large polarization reduction at high excita¬ 
tion powers. Fig. SJ^c) presents the spatial dependence of 
the electric field obtained from the approximate expres¬ 
sion Eq. m and it also gives values similar to those of the 
more general calculation. 

Note finally that, because the spin-spin coupling does 
not affect the spatially-averaged polarization < ^ >, 

the currents J _^ and _defined in Eq. [6] should also 

generate an increase of ^ at some value of r. This in¬ 
crease is not observed in the sample studied here because 
the magnitude of the spin-spin coupling is proportional 
to the polarization itself, which is small. A decrease of 
charge lifetime or an increase of spin lifetime will increase 
^ and should reveal an absolute maximum in the polar¬ 
ization at some distance from r = 0. As shown in the 
supplementary information, the maximum polarization 
may even exceed ^i. 

In conclusion, it has been shown both theoretically and 
experimentally that minority electron spin diffusion in 
the presence of slower diffusing photoholes generates a 
coupling between electon spins + and — such that the 
outward flow of spins ± generates an inward flow of spins 
T- This is a pure spin coupling which does not affect 
charge diffusion. The diffusion constant is then described 
by a matrix with negative nondiagonal elements and in¬ 
creased values of the diagonal elements with respect to 
the unipolar regime. This effect strongly reduces the spin 
polarization at the excitation spot for excitation powers 
that are sufficiently high to ensure that the inward am¬ 
bipolar spin currents are comparable with the outward 
diffusive currents. One of us (F. C.) is grateful to CON- 
ICYT Grant Becas Chile for supporting his work. 
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